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Abstract. Our understanding of the evolution of ^He and ''Li depends 
critically on the available data for these two elements at low metallicity. 
In particular, the degree to which there is a slope in an abundance vs 
metallicity regression can help determine the evolution of He, C, N, and 
O in dwarf galaxies in the case of ^He, and cosmic-ray induced nucleosyn- 
thesis of LiBeB in our own galaxy in the case of ^Li. Recent data and 
their implications will be discussed. 



1. Introduction 

There is a relatively large set of data available on the big bang nucleosynthe- 
sis (BBN) element isotopes of ^He and ^Li at low metallicity. It is common, 
particularly in the case of ^He, to perform a linear regression on the data with 
respect to some metallicity tracer such as 0/H or Fe/H. The intercept of such a 
regression can be directly related to the primordial abundance of that isotope, 
and the slope of the regression offers important clues to the nature of its chem- 
ical evolution. While one can not necessarily justify a linear relation from first 
principles, generally due to the quality of the data at low metallicity, such an 
approximation is acceptable. In fact, using the currently available ^He and ^Li 
datcj^, it is easy to show that a linear regression is significantly better than a 
weighted mean, yet more complicated fits using additional parameters generally 
do not yield a statistically significant improvement in the fit. 

Our inferences of the primordial abundances and evolution of the light el- 
ements are clearly tied to the quality of the data and our understanding of the 
systematic uncertainties in the derived abundances. Evolution is one the effects 
which is responsible for systematic uncertainties. In the case of ^He, the envi- 
ronment of the HII system is not pristine and includes non-primordial ^He. In 
addition, the true elemental abundances of ^He may be clouded due to effects 
such as underlying stellar absorption, collisional excitation, or flourecence. In 
the case of '^Li, abundances are contaminated by the non-primordial contribu- 
tion of ''Li from galactic cosmic-ray nucleosynthesis (GCRN) and uncertainties 
concerning the degree of stellar depletion of '^Li in pop II, halo stars. 



^Using, for example, the data of Izatov and Thuan (1998) on *He and Ryan, Norris, & Beers 
(1999) on ^Li. 
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2. ^He 

The ^He abundance has been determined from observations of Hell —>■ Hel re- 
combination hnes in a large sample of extragalactic HII regions (Pagel et al. 1992; 
Skillman & Kennicutt 1993; Skillman et al. 1994; Izotov, Thuan, & Lipovetsky 
1994,1997; Izotov k. Thuan 1998). Since ^He is produced in stars along with 
heavier elements such as Oxygen, it is expected that the primordial abundance 
of ^He can be determined from the intercept of the correlation between Y and 
0/H, namely Yp = y(0/H — > 0). A detailed analysis of the combined data 
(Olive, Skillman, k. Steigman 1997; Fields & Olive 1998) found an intercept 
corresponding to a primordial abundance 

Yp = 0.238 ± 0.002 ± 0.005 (1) 

The first uncertainty is purely statistical and the second uncertainty is an esti- 
mate of the systematic uncertainty in the primordial abundance determination. 
The helium abundance used in this analysis was determined using electron den- 
sities n obtained from SIX data. Izotov, Thuan, &; Lipovetsky (1994,1997) and 
Izotov & Thuan (1998) proposed a method based on several He emission lines to 
"self-consistcntly" determine the electron density. This method yields a higher 
primordial value 

Yp = 0.244 ± 0.002 ± 0.005 (2) 

Our interpretation of the evolution of ^He depends heavily on the slope 
of the ^He abundance with respect to a tracer element such as 0/H and/or 
N/H. While models of chemical evolution tend to give relatively low slopes 
(Ay/A(0/i/) ~ 20 — 60), the He data based on SII densities gives a much 
larger slope {AY/A{0/H) ~ 110 it 25), whereas the self-consistent method 
gives {AY/A{0/H) ~ 47 ± 26). The model calculations (Fields k Olive 1998, 
and references therein) depend crucially on the assumed yields of N in the AGB 
phase and on assumptions concerning hot-bottom burning. Many of the models 
attempting to reproduce the higher He slopes also rely on significant amounts 
of outflow in these dwarf galaxies. 

As can be ascertained from the brief discussion above, the method of anal- 
ysis has a huge impact on both the determination of the primordial ^He abun- 
dance and the slope of the He vs 0/H regression. Therefore, rather than discuss 
specific chemical evolution models in detail here, I will discuss some of the key 
sources of the uncertainties in the He abundance determinations and prospects 
for improvement. 

The He abundance is always quoted relative to H, e.g.. He line strengthes 
are measured relative to The H data must first be corrected for underlying 
absorption and reddening. Beginning with an observed line flux F{X), and an 
equivalent width W{X), we can parameterize the correction for underlying stellar 
absorption as 

- ^(^)(^) 

The parameter a is expected to be relatively insensitive to wavelength. A red- 
dening correction is applied to determine the intrinsic line intensity 1(A) relative 
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to Hp 



Xr{X) 



/(A) XAiX) 



(4) 



I {Hp) Xa{Hi3) 

where /(A) represents an assumed universal reddening law and C{H0) is the 
correction factor to be determined. By comparing Xii{\) to theoretical values, 
Xt{)^)^ wc determine the parameters a and C{HP) self consistently (Olive & 
Skillman, 2000), and run a Monte Carlo over the input data to test the robustness 
of the solution and to determine the systematic uncertainty associated with these 
corrections. 

In Figure 1 (from Olive &: Skillman 2000) , I show the result of such a Monte- 
Carlo based on synthetic data with an assumed correction of 2 A for underlying 
absorption and a value for C{H(3) = 0.1. The synthetic data were assumed 
to have an intrinsic 2% uncertainty. While the mean value of the Monte-Carlo 
results very accurately reproduces the input parameters, the spread in the values 
for a and C{Hf3) are generally a factor of 2 larger than one would have derived 
from the direct solution due to the covarience in a and C{Hf3). 
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Figure 1. A Monte Carlo determination of the underlying absorption 
a (in A ), and reddening parameter C{Hf5), based on synthetic data. 

The uncertainties found for Hf3 must next be propagated into the analysis 
for ^He, for which we follow an analogous procedure to that described above 
(Olive &: Skillman 2000). We again start with a set of observed quantities: line 
intensities /(A) which include the reddening correction previously determined 
along with its associated uncertainty which includes the uncertainties in C{H(3)] 
the equivalent width W{X); and temperature t. The Helium line intensities are 
scaled to Hfi and the singly ionized helium abundance is given by 



y+(A) 



/(A) E{HI3) ^{\) + a\ 



I{Hf3) E{\) ' W{\) ^(l + 7)/(r) 



(5) 



where E{X) / E{H (3) is the theoretical emissivity scaled to Hf3. The expression 

(5) also contains a correction factor for underlying stellar absorption, parame- 
terized now by a', a density dependent collisional correction factor, (1 -|- 7)"^, 
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and a flouiccence correction which depends on the optical depth r. Thus y"*" 
imphcitly depends on 3 unknowns, the electron density, n, a', and r. 

One can use 3-6 lines to determine the weighted average helium abundance, 
y. From y, we can calculate the deviation from the average, and minimize 
X^, to determine n,a', and r. Uncertainties in the output parameters are also 
determined. 

This procedure differs somewhat from that proposed by ITL, in that the 
above is based on a straight weighted average, where as ITL minimize the 
difference of a ratio of He abundances (to one wavelength, say A4471) to the 
theoretical ratio. When the reference line is particularly sensitive to a systematic 
effect such as underlying stellar absorption, this uncertainty propagates to all 
lines this way. In our case, the individual uncertainties in the line strengths are 
kept separate. 

Finally, as in the case for the hydrogen lines, we have performed a Monte- 
Carlo simulation of the data to test the robustness of the solution for n, a', and 
r (Olive Sz Skillman 2000). In Figure 2, I show the result of a single case based 
on the data of Izatov and Thuan (1998) for SBS1159-I-545. Here, the helium 
abundance and density solutions are displayed. The vertical and horizontal lines 
show the position of the IT solution. The circle shows the position of the our 
solution to the minimization, and the square shows the position of the mean 
of the Monte-Carlo distribution. The spread shown here is significantly greater 
than the uncertainty quoted by IT. 
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Figure 2. A Monte Carlo determination of the helium abundance 
and electron density (in cm~^). Solutions for a' and r are not shown 
here. 



3. ^Li 

The population II abundance of ^Li has been determined by observations of over 
100 hot, halo stars, and is found to have a very nearly uniform abundance (Spite 
& Spite 1982). For stars with a surface temperature T > 5500 K and a metal- 
licity less than about l/20th solar, the abundances show little or no dispersion 
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beyond that which is consistent with the errors of individual measurements. The 
Li data from Bonifacio Sz Molaro (1997) indicate a mean ^Li abundance of 



The small error is statistical and is due to the large number of stars in which 
'''Li has been observed. 

There is, however, an important source of systematic error due to the pos- 
sibility that Li has been depleted in these stars, though the lack of dispersion 
in the Li data limits the amount of depletion. In fact, as discussed by Scan 
Ryan (these proceedings, and Ryan, Norris, &; Beers, 1999, hereafter RNB) a 
small observed slope in Li vs Fe and the tiny dispersion about that correlation 
indicates that depletion is negligible in these stars. Furthermore, the slope may 
indicate a lower abundance of Li than that in (6) . 

For reference, the weighted mean of the ^Li abundance in the RNB sample 
is [Li] = 2.12 ([Li] = log '^Li/H + 12). It is common to test for the presence of a 
slope in the Li data by fitting a regression of the form [Li] = a + P [Fe/H]. The 
RNB data indicate a rather large slope, (3 = 0.07 — 0.16 and hence a downward 

shift in the "primordial" lithium abundance A [Li] = —0.20 0.09. Models of 

galactic evolution which predict a small slope for [Li] vs. [Fe/H], can produce a 
value for (3 in the range 0.04 - 0.07 (Ryan et al. 2000). Of course, if we would 
like to extract the primordial ^Li abundance, we must examine the linear (rather 
than log) regressions. For Li/H = a' + 6'Fc/Fcq, wc find a' = 1 — 1.2 x 10^^*^ 
and 6' = 40 — 120 x lO"^'^. A similar result is found fitting Li vs O. Overall, 
when the regression based on the data and other systematic effects are taken 
into account a best value for Li/H was found to be (Ryan et al. 2000) 



1.23 X 10"^'^. The linear slope produced by the model is 6' = 65 x 10"^'', 
and is independent of the input primordial value (unlike the log slope given 
above). The model (discussed in detail in Fields & Olive, 1999a, b and Fields 
et al. 2000) includes in addition to primordial ^Li, lithium produced in galactic 
cosmic ray nucleosynthesis, (primarily a + a fusion) in addition to 'Li produced 
by the z^-process during type II supernovac. As one can see, these processes are 
not sufficient to reproduce the population I abundance of ''Li, and additional 
production sources are needed (see e.g. Matteucchi, these proceedings). 

4. Concordance 

Bearing in mind the degree of uncertainty in the derived primordial abundances, 
one can test the concordance of ^He and ^Li with the prediction of BBN. This 
is best summarized in a comparison of likelihood functions as a function of the 
one free parameter of BBN, namely the baryon-to-photon ratio 77. By combin- 
ing the theoretical predictions (and its uncertainties) with the observationally 
determined abundances discussed above, we can produce individual likelihood 
functions (Fields et al. 1996) which are shown in Figure 4. A range of primordial 



Li/H = (1.6 ±0.1) X 10- 



(6) 




(7) 
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Figure 3. Contributions to the total predicted lithium abundance 
from the adopted GCE model of Fields & Olive ( 1999a, b), compared 
with low metallicity stars (RNB) and a sample of high metallicity stars. 
The solid curve is the sum of all components. 

^Li values are chosen based on the analysis in Ryan et al. (2000). The double 
peaked nature of the ^Li likelihood functions is due to the presence of a mini- 
mum in the predicted lithium abundance in the expected range for r/. For a given 
observed value of ^Li, there are two likely values of 77. As the lithium abundance 
is lowered, one tends toward the minimum of the BBN prediction, and the two 
peaks merge. Also shown are both values of the primordial ^He abundances dis- 
cussed above. As one can see, at this level there is clearly concordance between 
"^He, ^Li and BBN. 



5. LiBeB 

The production of ^Li by galactic cosmic-ray nucleosynthesis shown in Figure 3, 
is accompanied by the production of the heavier intermediate elements Be and B 
(Reeves, Fowler, &: Hoyle 1970; Meneguzzi, Audouze, & Reeves 1971). Standard 
GCRN is dominated by interactions originating from accelerated protons and 
a's on CNO in the ISM, and predicts that BeB should be "secondary" versus the 
spallation targets, giving Be a (Vangioni-Flam, Casse, Audouze, k. Oberto 
1990). However, this simple model was challenged by the observations of BeB 
abundances in Pop II stars, and particularly the BeB trends versus metallicity. 
Measurements showed that both Be and B vary roughly linearly with Fe, a so- 
called "primary" scaling. If O and Fe are co-produced (i.e., if 0/Fe is constant) 
then the data clearly contradicts the canonical theory, i.e. BeB production via 
standard GCR's. 

These observations led to the creation of many new models of cosmic-ray 
nucleosynthesis (Casse, M., Lehoucq, R., & Vangioni-Flam, E. 1995; Vangioni- 
Flam et al. 1996; Ramaty et al. 1997, 1999) which include a dominant primary 
component of BeB. Such models were discussed here by Casse, Parizot and 
Ramaty. 
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Figure 4. Likelihood distributions for four values of primordial ^Li/H 
(lO^^x ■^Li = 1.9 (dashed), 1.6 [dotted), 1.23 (solid), and 0.9 {dash- 
dotted)), and for ^He (shaded) for which we adopt Yp = 0.238 it 0.002 it 
0.005 (Eq. (1)). Also shown by the long dashed curve is the likelihood 
function based on the "^He abundance from Eq. (2). 



As was discussed here by Deliyannis and Israelian, there is growing evidence 
that the 0/Fe ratio is not constant at low metallicity (Israelian, Garcfa-Lopez, 
& Rebolo 1998; Boesgaard et al. 1999a), but rather increases towards low metal- 
licity. This trend offers another solution to resolve discrepancy between the ob- 
served BeB abundances as a function of metallicity and the predicted secondary 
trend of GCR spallation (Fields & Olive 1998). As noted above, standard OCR 
nucleosynthesis predicts Be cx O^, while observations show Be ~ Fe, roughly; 
these two trends can be consistent if 0/Fe is not constant in Pop II. A com- 
bination of standard GCR nucleosynthesis, and i/-process production of ^^B is 
consistent with current data. 

The determination of abundances from raw stellar spectra requires stel- 
lar atmosphere models. The atmospheric models require key input parameters, 
notably the effective temperature Tefj and surface gravity g, and assumptions 
regarding the applicability of local thermodynamic equilibrium (LTE). Unfortu- 
nately, there is no standard set of stellar parameters for the halo stars of interest. 
In practice, different groups derive abundances via different procedures, which 
give similar results but retain systematic differences. The systematic differ- 
ences in the data can in fact obscure the BeB-OFe trends one seeks. Thus, to 
derive meaningful BeB fits, one must systematically and consistently present 
abundances derived under the same assumptions and parameters for stellar at- 
mospheres. 

It is not possible to overly stress the importance of reliable stellar data. 
The Balmer-line data appear to be self-consistent, and are probably the most 
reliable. However, because other scales such as those based on the IRFM are 
commonly used, I would like to point out that there are significant differences in 
the reported data. To illustrate the point consider for example the case of the 
star BD 3° 740. From Axer et al. (1994), whose data is based on the Balmer 
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line method we find this star to have {T^sMg, [Fe/H]) = (6264, 3.72, -2.36). 
The beryllium and oxygen abundances for this star was reported by Boesgaard 
et al. (1999b) and (1999a). When adjusted for these stellar parameters, we 
find [Be/H] = -13.36, and [0/H] = -1.74. In contrast, the stellar parameters 
from Alonso et al. (1996a) based on the IRFM (IRFMl) are (6110,3.73,-2.01) 
with corresponding Be and O abundances of -13.44 and -2.05. Garcia-Lopez 
et al. 1998 use a calibrated IRFM (IRFM2) based on Alonso et al. (1996b) 
and take (6295,4.00, -3.00). For these choices, we have [Bc/H] = -13.24 and 
[0/H] = -1.90. Notice the extremely large range in assumed mctallicities and 
the difference in the two so-called IRFM temperatures. While this may not be 
a typical example of the difference in stellar parameters, it is differences such 
as this (and this star is not unique) that accounts for the difference in our 
results and the implications we must draw from them. Uncertainties in [Fe/H] 
in particular, make modeling extremely difficult. This is especially true of one 
attempts to model the correlations of BeBO with respect to Fe/H. 

Below, we will present results based for the available BeBOFe data based 
on three methods of analysis. Wc will refer to these as the Balmcr line data and 
the IRFMl, 2 data. Complete results of this analysis can be found in Fields et 
al. (2000). There are a total of 36 stars with low metallicity OH data. Of these, 
roughly 2/3 have available data using one of the systematic methods described 
(Balmer, IRFMl, IRFM2). In each case, one finds a significant slope for [0/Fe] 
vs [Fe/H] ranging from -0.32 to -0.51. 

Of key importance to the modeling of the BeB evolution is the determination 
of a primary or secondary source for the BeB isotopes. Primary vs secondary is 
typically ascertained by fitting the BeB data versus a tracer element. Histori- 
cally, Fe/H was used even though the actual production of LiBeB is independent 
of [Fe/H]. This is justifiable so long as [0/Fe] is constant. As one can see in 
the tables below, the data seem to indicate that Be is mostly primary with re- 
spect to Fe, and secondary with respect to 0/H. (IRFM2 should be considered 
suspect as the derived parameters were obtained outside the limits of validity 
of the calibration.) This is what one would expect if [0/Fe] is not constant as 
the OH data now indicate. B on the other hand shows primary evolution with 
respect to [0/H] and an even flatter evolution with respect to Fe/H. This too is 
expected if the i^-process plays a significant role in the production of ^^B. 



Table 1. Slopes for Be versus Fe and O. 



method 


number 


tracer 


slope 


number 


tracer 


slope 


Balmer 


22 


Fe 


1.39 ±0.16 


19 





1.78 ±0.19 


IRFMl 


22 


Fe 


1.23 ±0.14 


21 





1.83 ±0.19 


IRFM2 


18 


Fe 


1.18 ±0.11 


18 





1.36 ±0.09 



The models for primary and secondary production of BeB are all physical. 

What is unclear however, is which is dominant over the history of the Galaxy 
and at what epoch. If both mechanisms are operative, it is reasonably certain 
that primary mechanisms should dominate in the early Galaxy and that sec- 
ondary mechanisms should dominate later. The cross-over or break point can 
be determined from the data (in principle) by fitting to both linear and quadratic 
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Table 2. Slopes for B versus Fe and O. 



method 


number 


tracer 


slope 


number 


tracer 


slope 


Balmer 


11 


Fe 


0.78 ±0.22 


10 





1.23 ±0.32 


IRFMl 


9 


Fe 


0.73 ± 0.19 


9 





0.98 ± 0.28 


IRFM2 


11 


Fe 


0.72 ± 0.14 


11 





1.02 ±0.16 



components, 



A 
H 



o 



ai 



0/H 



+ a2 



( 0/H 



(O/H)0 ^V(0/H)e 



(8) 



for A € BcB. The resulting coefficients and break points for Be are found in 
Table 3. As one can see, for Balmer and IRFMl, the break point occurs at 
low [0/H], indicating that most of the evolution in the observed data has been 
secondary. To fully resolve this issue, a larger and systematic data set is required. 



Table 3. Break points for Be versus O. 



number 


method 


a\ 


Oi2 


[O/HJeq 


19 


Balmer 


0.042 ± 0.003 


2.30 ± 0.70 


-1.75 


21 


IRFMl 


0.034 ± 0.034 


2.14 ±0.53 


-1.79 


18 


IRFM2 


0.111 ±0.031 


2.57 ± 0.76 


-1.37 



Finally, in Figure 5 (from Fields et al. 2000), the evolution of BeB with 
respect to 0/H is shown in a simple closed box model of chemical evolution. 
In addition to standard GCR nucleosynthesis, a primary component based on 
the superbubble accelerated particle spectrum of Bykov (1999) is included along 
with the neutrino-process for ^^B (Fields et al. 2000). The secondary GCR 
cosmic-ray flux is normalized by the solar abundance of Be and is consisitent 
with the present cosmic-ray flux scaled by the star formation rate. The break 
point (from Table 3) determined the relative scaling of the primary component 
to the secondary one, and the neutrino process is scaled to the solar ^^B/^°B 
ratio. 
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Figure 5. Be vs O {top panel) and B O {bottom panel). Data 
shown are the Balmer points, which are found to have a break point 
as indicated. Models are adjusted to have the break point and 0/Fe 
slope of these data. 
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